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ABSTRACT 

We report the detection of hydrogen chloride, HC1, in diffuse molecular clouds 
on the line of sight towards the star-forming region W31C (G10.6-0.4). The 
J = 1 —0 lines of the two stable HC1 isotopologues, H 35 C1 and H 37 C1, are observed 
using the lb receiver of the Heterodyne Instrument for the Far-Infrared (HIFI) 
aboard the Herschel Space Observatory. The HC1 line is detected in absorption, 
over a wide range of velocities associated with diffuse clouds along the line of 
sight to W31C. The analysis of the absorption strength yields a total HC1 column 
density of few 10 13 cm -2 , implying that HC1 accounts for ~ 0.6 % of the total gas- 
phase chlorine, which exceeds by a factor of ~ 6 the theoretical model predictions. 
This result is comparable to those obtained from the chemically-related species 
H2CI 4 " and HC1 + , for which large column densities have also been reported on the 
same line of sight. The source of discrepancy between models and observations 
is still unknown; however, the detection of these Cl-bearing molecules, provides 
key constraints for the chlorine chemistry in the diffuse gas. 

Subject headings: astrochemistry — submillimeter: ISM - ISM: molecules- 
ISM:abundances 
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1. Introduction 

Hydride molecules play an important role in interstellar chemistry, as they are often 
stable end points of chemical reactions, or represent important intermediate stages of the 
reaction chains theorized to form gas-phase molecules. This makes hydrides a sensitive test 
of these chemical models, as well as potential tracers of other molecules of interest e.g. 
molecular hydrogen. Due to their small moment of inertia, hydrides have their fundamental 
rotational lines in the submillimeter band. The Heterodyne Instrument for the Far-Infrared 
(HIFI) on board the Herschel Space Observatory is providing invaluable data on interstellar 
chemistry in general, and in particular on hydride molecules within the Milky Way and local 
galaxies. Halogen atoms fluorine (F) and chlorine (CI), with estimated solar abundances of 
3.6 x 10~ 8 and 3.2 x 10~ 7 relative to hydrogen (Asplund et al. 2009), are of special interest 
because they are the major atoms (neutral F and ionic Cl + ) in diffuse environments that 
react exothermically with molecular hydrogen to form hydride molecules (HF and HC1). 
HF and HC1 are strong bound systems, only destroyed by photodissociation, reactions with 



He + , H^~, and C + and, in the case of HC1, also by photoionizatio n. Chlorin e 



been determined by extensiye theoretical and observational work ( 



1974 



van Dishoeck fc Black 



1986 



Blake et al. 



1986 



Schilke et al. 



ura 



1974; 



1995 



chemistry has 



Dalearno et al. 



Federman et al. 



19951 ). In diffuse c louds, CI a toms can be ionized by UV photons at wavelengths between 



91.2 and 95.6 nm ( jJura 



19741 ) . The resulting ion, Cl + , reacts rapidly exothermically (by 



0.17 eV) with H 2 to form HC1+ 



cr 



h 2 ->■ Hcr 



H. 



which in turn reacts with H 2 to form chloronium, H 2 C1 



HC1+ + H, -> H 2 C1+ + H. 



(2) 
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HC1 is then formed through dissociative recombination (DR) of H2C1 + 



H 2 C1 + + e -> 




(3) 



From reactions (1) - (3), HC1 is expected to be abundant in regions containing both 
H2 and chlorine-ionizing photons. 

Observations of HC1 have been limited to ground-based observatories at high-altitude 
sites, under extreme good weather conditions, or space-missions (due to its large 
rotational constant which places the ground state transition, at 625.9187 GHz, near a 
strong atmospheric water absorption line). The first detections of interstellar chlorine- 
containing molecu les were obtained using the NASA's Kuiper Airborne Observatory (KAO). 



Blake et al. 



(119851 ) observed the HC1 ground-st ate J = 1 — rotat i onal t ransition in emission 



in OMC-1 with angular resolution of about 2'. 



Zmuidzinas et al. 



( 19951 ) detected the line in 



absorption towards Sagittarius B2. Following these results, ground-based observations have 



lead to the detection of the H 35 C1 and 



stars and active star-forming regions (jSchilke et al 



2010|). 



H 37 C1 isotopo 


ogues 


in environments such as evo 


(Schilke et al. 


1995 




Salez et al. 


1996 




Pene et al. 



ved 



CI chemistry in diffuse clouds on the other hand, has not been fully constrained 
because of the paucity of observations of Cl-bearing species in the interstellar medium. 
Prior to Herschel/HIFI, HC1 in diffuse clouds was observed mainly by m eans of ultraviole t 



absor ption studies leading only to a tentative detection towards ( Oph ( IFederman et al. 



19951 ). with an HC1 column density of 2.7 x 10 11 cm 2 . Given an atom ic chlorine column 



1995|), the jV(HCl)/JV(Cl) 



density of 3.0 x 10 14 cm 2 for this line of sight (IFederman et al. 
ratio is 9 x 10~ 4 . This result is in goo d agreement with the theoretical models o f the CI 



chemistry (Ivan Dishoeck &: Black 



19861 and more recently 



Neufeld fc Wolfire 



20091 . hereafter 

NW09). The NW09 models predict column densities in the range of 10 10 - 10 12 cm -2 for 
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HC1, H2C1 + and HC1 + , for a diffuse molecular cloud of density nn = 10 2 5 cm 3 and xuv 



in the range of 1 - 10, where \ 



mean interstellar value f lDraine 



JV 



is the UV radiation field normalized with respect to the 



1978 



see Figures 6 and 7 in NW09). Model predictions also 



identified H2CI 4 " and HC1 + as potentially detectable Cl-bearing species. 

The Herschel/HIFI instrument has indeed allowed observations for the first time of 
these two new Cl-bear ing species in the interstellar medium , as well as the detection of 



HC1 in diffuse cloud s ( iLis et al 



2010; 



De Luca et al. 



( iCodella et al. 



20121 ) and in protostellar shocks 



20121 ) . The molecular ions H2C1 + and HC1 + have been detected in diffuse 



clouds on the lines of sight to bright submillimeter continuum s ources. H?C1 + has been 



detected in absorption towards NGC 63341, Sagittarius B2 (S) (|Lis et al. 



Sgr A and in emission towards Orion Bar, Orion South ( jNeufeld et al. 



20101 ). W31C, 



2OI2I). The H 2 C1+ 



column densities obtained from these studies, imply that chloronium accounts for ~ 

4 - 12 % of chlorine nuclei in the gas phase. This result is at least a factor ~ 10 larger than 

that predicted by the chemical models, which pr edict a H9CI 4 " and H C1 + joint contribution 



of ~ 2 % to the total gas-phase chlorine budget (jNeufeld et al. 



2012 



). Similar resu 



been derived from observations of HC1 + towards W31C and W49N ( De Luca et al. 



■;s have 



2OI2J) 



where HC1 + has been detected in absorption with large column densities suggesting a 
3 - 5 % contribution of this species to the total gas-phase chlorine content. 

The discrepancy between the chlorine chemistry models and observations is very 
puzzling; diffuse clouds contain only simple molecules, and the number of reactions involved 
in describing their abundances is fewer than in dense molecular clouds. Therefore, a test of 
the basic interstellar chemistry should be possible, establishing much higher standards for 
the modeling of diffuse clouds compared to those of dense clouds. The diffuse cloud models 
thus should be able to reproduce the measured abundances to a factor of two or better 



(Ivan Dishoeck 



19901 ) 
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In this paper, we report the result of Herschel/HIFI observations of HC1 J = 1 — on 
the line of sight towa rds W31C (G10.6-0.4). These observations were previously attempted 



by 



Peng et al. 



(120 1Q[ ) using the Caltech Submillimeter Observatory (CSO), but yielded no 
HC1 detection due to the high rms of ~ 0.15 K at a velocity resolution of 0.31 km s" 1 in 
their observations. W31C is an extremely active region of high- mass star formation, and 
one of the three bright HII regions of the W31 complex, harb oring an extremely luminous 



submillimeter and infrared continuu m source ( 
a kinematic distance of 4.8lgg kpc fjFish et al. 



'TR 



10 7 L P 



0) 



Wright et al. 



19771 ). With 



20031 ). the W31C line of sight intersects 



several foreground molecular clouds from the Milky Way spiral arms. The source location 
and its strong continuum flux, make W31C one of the best sources towards which to carry 
out absorption line studies. We compare the Herschel/HIFI results with current chlorine 
chemistry models and with observations of other Cl-bearing molecules detected along the 
same line of sight and toward the galactic center source Sgr B2(S). 



Observations 



The observations were performed with the Herschel/HIFI instrument (jde Graauw et al. 



20101 ) on March 2 2010, as part of the PRISMAS (Probing Interstellar Molecules with 
Absorption lines Studies) guaranteed time key program. Both H 35 C1 and H 37 C1 J = 1 — 
lines, with rest frequencies of 625.918 and 624.977 GHz respectively, were observed 
simultaneously using the HIFI band lb receiver, with 3 shifted local oscillator (LO) settings 
and the wide band spectrometer (WBS). The WBS provides a spectral resolution of 1.1 
MHz over a 4 GHz Intermediate Frequency (IF) bandwidth. The dual beam switch (DBS) 
observing mode was used with reference beams located 3' on either side of the source 
position (a J20 oo = 18 h 10 m 28.700 s and 5 J20 oo = -19°55'50.00") along an East- West axis. 

The data have been reduced using the Herschel Interactive Processing Environment 
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(HIPE) flOtt et al 



2010() with pipeline version 5.2. The resulting Level 2 double sideband 
spectra were exported into FITS format for subsequent data reduction and analysis using 
the IRAM GILDAS packagd_|. The spectra obtained at different LO setting and for each 
polarization were inspected for possible contamination by emission lines from the image 
sideband and averaged to produce the final spectra shown in Figure [TJ Beam measurements, 
reported on November 17 of 2010, towards Mars at 610 GHz give a main beam (t]^) of 
0.744 and 0.764, for the Horizontal (H) and vertical (V) polarization, respectively. The full 
width at half maximum (FWHM) HIFI beam size at the HC1 J = 1 — frequency is ~ 33". 



3. Results 

The spectra obtained for the two stable hydrogen chloride isotopologues, H 35 C1 and 
H 37 C1 towards W31C are shown in Figure HJ upper and lower panels respectively. The data 
quality is excellent with a double sideband continuum antenna temperature of ~ 1.6 K, 
and a rms noise of ~ 6 mK. Figure [1] shows the H 35 C1 and H 37 C1 emission lines from the 
background source and a clear absorption feature in the H 35 C1 spectrum, at velocities 
~ 10 - 50 km s" 1 . The emission lines from the source will be discussed in an extended 
paper, which will include the HC1 observations towards all PRISMAS sources (Monje et al. 
in prep.). The H 35 C1 absorption feature is analogous to the one seen in t he spectra of other 



Neufeld et al 



20121 ) and corresponds 



molecules such as HF, CH and H 2 (see Figure 2 in 
to the foreground clouds on the line of sight towards W31C. The foreground clouds are 
known to harbor low density and cool molecular gas (Tki n ~ 50 - 70 K from H 2 and 
~ 100 K from H I), where molecular excitation can be highly subthermal, w ith excitation 



temperatures close to the cosmic microwave background radiation ~ 2.73 K (IGodard et al. 



x http: / /iram.fr/IRAMFR/GILDAS / 
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2 1 Oh . Consequently, molecular emission lines are very weak and these clouds are best 
studied through absorption spectroscopy 

The corresponding H 37 C1 absorption line is contaminated by interfering emission of 
dimethyl ether (CH3OCH3), methyl cyanide (CH 3 CN, v =0) and sulphur dioxide (SO2). Th e 



emiss ion lines are modeled with the program XCLASS (jSchilke et al. 



1998 



Comito et al. 



20051 ). which assumes local thermodynamic equilibrium (LTE) and uses th e rest frequen cies 



provided by the Cologn e Database for Molecular Spectroscopy (CDMS) (jMiiller et al 



2001 



Mill 



er et al 



2005 



) and the Jet Propulsion Laboratory (JPL) molecular spectroscopy 



database (IPickett et al. 



19981 ). The result from the line fits exhibits two relatively strong 
CH 3 CN and CH3OCH3 lines and a weak SO2 line within the H 37 C1 absorption spectrum, 
see Figure [2j To obtain the best LTE fit for CH3OCH3 and SO2 additional lines within 
the same LO tuning are used, see Figure [2j For the CH3CN additional transition lines 
at 533.123 GHz (J = 29 - 28) and 606.533 GHz (J = 33 - 32) are used, see Figure [3j 
The parameters used in the LT E model, shown in T able 1, assume that the CH 3 CN and 



SO2 are tracers of the hot core ( iBeltran et al. 



20 111 ) while the CH3OCH3 arises from the 



more extended envelope. However, the W31C region is known to be a complex source with 
outflow activity from embedded high-mass protostars, methanol masers, a dense rotating 
toroid and infalling material from the molecular envelope. The LTE emission line fit model 
is thus an approximation and tighter constraints using the combination of key molecules 
such as methanol, formaldehyde, CO, CS and SO and more complex chemical modeling 
are needed. Due to the uncertainties of the modeling, we only use the contamination-free 
spectra (grey shaded area in Figure [2]) to derive an estimate of the 35 C1/ 37 C1 isotopic ratio, 
as described below. 

The rotational levels of HC1 have a hyperfine splitting, caused by the interaction 
of the quadrupole moment of the chlorine nucleus (I = 3/2) and the electric field. The 
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J = — 1 transition is thus split into three components with relative strengths of 50 %, 
16.7 %, and 33.3 %, and with velocity offsets (Av h ) of 0, -6.3, and +8.2 km s" 1 for 
H 35 C1, and 0, -5 and 6.5 km s _1 for H 37 C1. The Avh is larger than the velocity difference 
between consecutive velocity components, producing overlaps between different velocity 
components that results in a complex profile. In order to determine the velocity structure 
of the foreground absorbing gas, we extract the signal associated with each HFS component 



Godard et al. 



(120 121 ) . In Figure HI we present 



using the numerical procedure described in 
the H 35 C1 spectrum, with flux normalized with respect to the single sideband continuum 
(2Tl/Tc - 1, where Tl/Tc is the line-to-continuum ratio) in the velocity range from vlsr 
~ 5 to 70 km s -1 . The decomposition of the 625 GHz line into three hyperfine components, 
the main F = 5/2 - 3/2 and the two satellite F = 1/2 — 3/2 and F = 3/2 — 3/2 hyperfine 
components, is also plotted in Figure SJ We add the three hyperfine components, shifting 
(in velocity) the satellite components accordingly to the main hyperfine component velocity 
shift, and obtain a hyperfine deconvolved spectrum, equivalent to a splitting-free spectrum, 
where each HFS component contributes to the optical depth. In order to estimate the line 
width and the center velocity of each individual velocity component, we fit the spectrum 
with a set of Gaussian curves (see Figure [5]). The five Gaussian components used in our fit 
were set to give a better correspondence to o ther molecular trans itions obse rved toward s 



the s ame background source, such as H 2 C1 + (INeufeld et al. 



20121 ) and CH (IGerin et al. 



20101 ). For comparison, we also plot in 



absor ption profiles from 



Neufeld et al 



'igure|5lthe H 2 35 C1+ and H 2 37 C1+ HFS deconvolved 



( 120121 ) and the H I absorption spectrum (IFish et al. 



20031 ). The H2 35 C1 + and H2 37 C1 + profiles are in good agreement with that of H 35 C1 for most 



velocity intervals, as expected from precursor species of HC1. 
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3.1. Column densities 



We calculate the hydrogen chloride column densities for the velocities intervals 
associated with the different Gaussian components. First, we derive optical depths of the 
HC1 lines (r = -ln[2T^/Tc-l]), assuming that the foreground absorption completely covers 



the continuum source and that all HC1 molecules are in the ground state . We derive t 



lie 



Neufeld et al. 



H 35 C1 column densities for each LSR velocity range using equation 3 of 
(120101 ). where the absorption optical depth for the H 35 C1 J = — 1 transition, integrated 
over velocity, is given by 

A u ig u \ 3 



rdv 



-iV(H 35 Cl) => iV(H 35 Cl) = 6.5211 x lO 1 ^ rdv [cm' 2 ] 



where g u — 3 and gi = 1 are the degeneracies of the upper (J = 1) and lower (J = 0) 
state. Since we use the deconvolved spectrum without hyperfme splitting, A = 478.96 /xm 
is the transition wavelength, and A u \ = 1.17 x 10 -3 s~ x is the spontaneous radiative decay 
rate. Table 2 shows the HC1 column densities and the abundances with respect to the total 
hydrogen column density. We calculate the 35 CI/ 3 CI ratio in the velocity interval free from 
emission line contamination (shaded interval in Figure |2]), which corresponds to velocities 
from 38 to 40.5 km s _1 . Th e 35 C1/ 37 C1 ratio obtained is equal to 2.9, in good agreement 



19891 ). Thus, the HC1 abundances are 



with the solar value of 3.1 (lAnders fc Grevessd 
calculated with respect to the total molecular and neutral hydrogen assuming a 35 C1/ 37 C1 
of 3.1 in all velo c ity in tervals. The H I column densities in the foreground ga s are obtaine d 



from 



Fish et al. 



20101) 



(120031 ). while the molecular hydrogen is obtained from CH (jGerin et al. 



assuming a 



10- 8 (ISheffer et al. 



inear correspondence between CH and H2 with a scaling factor of 3.5 x 



20081 ). The average HC1 fractional abundance with respect to the total 
hydrogen is 6.1 ± 0.6 x 10~ 10 . Using a chlorine abundanc e in di ffuse clouds, N(Cl)/Nn, 



of 1.03 x 10 7 measured by 



Moomey. Federman fc Shefferl (120121 ) towards nearby (a few 



hundred pc) stars, we obtain that the average fraction of gas-phase CI in HC1 is ~ 0.6 %. 
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The uncertainties in the column densities are the random noise and the systematic errors 
introduced by the calibration uncertainties and the robustness of the HFS decomposition 
procedure. The 1-er uncertainty in the column density is dominated by the error in the 
linear baseline fitting. The total HIFI calibration uncertainties for band lb are < 9.3 %, 
which includes the contribution from the sideband gain ratio, beam efficiency, pointing, hot 
and cold coupling, where all errors ar e added in quadr a ture. Detailed information about 



the HIFI calibration can be found on 



Roelfsema et al 



(12012). The error intr oduced by 



the robustness of the deconvolution procedure is ~ 14 % (IGodard et al. 



20121 ). Adding in 



quadrature all the errors, the total uncertainty in the column densities is about 19 



Chemical model 



We use the Meudon PDR code ( ILe Petit et al 



20061 ) to compute the fractional 



abundances of the various chlorine bearing species, i.e. the abundance of each species 
normalized relative to the chlorine elemental abundance, as a function of the total extinction 
of the cloud, A Vyto t- The Meudon PDR code is a one- or two-sided model, where the 
molecular cloud is modeled as a stationary plane-parallel slab of gas and dust, exposed 
to incident radiation field from a bright star or the Interstell ar Standard Radiatio n Field 



Neufeld et al 



(120121 ). which 



(ISRF). As input parameters, we use the best fit model from 
corresponds to an initial density of Tin = 316 cm -3 , a radiation field intensity from both 
sides of 10 and a primary ionization rate of 1.810 x 10 ~ 16 s~ l . Using the standard Galactic 



20031 ). A V)to t probed is 



gas-to-dust ratio Nn/A v = 1.93 x 10 21 mag" 1 cm -2 (IWhittet 
~ 3 - 5. The obtained abundances normalized relative to the relevant gas-phase elemental 
abundance are shown in Figure 6-left, the results indicate that the HC1 + and H 2 C1 + account 
jointly at best for ~ 2 % of the elemental chlorine in the gas phase while HC1 accounts for 



0.1 %. The modeled HC1 fractional abundance under predicts the results obtained from 
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the observations by an average value of about 6. 



The DR studies of H 2 C1 + and the corresponding product branching ratios are not yet 
available from laboratory experiments. In their models, NW09 assumed a branching ratio 
of 90 % for CI and 10 % for HC1, motivated by the low HC1 abundance observed towards 
the diffuse cloud £ Oph. We have thus investigated the effect of the branching ratio of 
H2C1 + DR for the physical conditions described above, based on the correlation obtained 
between the branching ratio of triatomic dihydride ions fo r 3-bo dy dissociation and the 
corresponding energy r elease obtained by 



Roueff &: Herbst 



Roueff fc Herbst 



(120111 ). From that correlation, 



(120111 ) obtain a branching ratio of 56 % for the 3-body dissociation of 
H 2 C1 + . We estimate the fractional abundances corresponding to the new branching ratio 
and plot them as dashed lines in Figure 6-right, in comparison with those obtain with 
previous assumption (90 % for CI) in full lines, as a function of A Vtot . The results show that 
differences occur only for the minor species, i.e. HC1, HC1 + and H 2 C1 + while the reservoirs 
CI and Cl + are not affected. The increase in the HC1 column density is about a factor of 
~ 5, in closer agreement with our observational results. However, the increase for HC1 + 



and H 2 C1 + is much smaller, and then this result cannot explain all the discrepancies. 



5. Discussion 



We present the detection of H 35 C1 and H 37 C1 J = — 1 transition in absorption on 
the line of sight towards the bright submillimeter source W31C. The obtained HC1 column 



density 



B2 (S ) ([Lis 



2 x 10 cm ), is comparable with that obtained towards the lin e of sight to Sgr 



2012|) and 



et al 



20 1 01) . Given the HC1 + and H 2 C1 + column densities from 



De Luca et al. 



Neufeld et al. 



(J20JJ), the HC1+/H 2 C1+ and HC1/H 2 C1+ column density ratios, of 
~ 1 and ~ 0.3 respectively, are within the range predicted by the models of diffuse clouds. 
However, the averaged HC1 abundance with respect to the total hydrogen content of 6.1 ± 



13 



0.6 x 10~ 10 suggests a HC1 fractional abundance with respect to the chlorine elemental 
abundance (N(CI)/Nk = 1.03 x 10~ 7 ) a factor of about 6 larger than those predicted by 
chemical models of diffuse clouds. Similar discrepancies between models and observations 
were also found for other Cl-bearing gas-phase spe cies, HC1 + and H 2 C1 + towards W31 C, 



Sgr A (+50 km s 1 molecular cloud) and W49N (ILis et al. 



2010 



De Luca et al. 



2012 



Neufeld et al. 



20121 ). 



Lis et al. 



(120 lOf ) suggested a geometrical explanation for the discrepancy between the 
observed column densities and the models, where enhancement in the absorbing column 
density could result from multiple PDRs present along the line of sight. This scenario could 
explain higher column densities, but the HC1/C1 ratio will remain constant and hence this 
solution is not likely to reproduce the W31C line of sight results. 



Neufeld et al 



( 120121 ) explored several possible explanations for the deficiencies in the 
model predictions, e.g. the rate coefficient assumed for dissociative recombination (DR) 
of HC1 + and uncertainties in the assumed dipole moment of H2C1 + , none of which proved 
satisfactory. We have investigated the effect of varying the branching ratio of H 2 C1 + DR. 
By taking a branching ratio of 56 % for the chlorine (instead of 90 % assumed in NW09 
models), the predicted fractional abundance of HC1 with respect to the gas-phase chlorine 
abundance increases by about a factor of ~ 5. However, laboratory results are needed 
to confirm the branching ratio and validate this argument, along with the study of other 
unknowns in the chlorine chemistry, such as the DR rates, and the dipole moment of H 2 C1 + . 



HIFI has been designed and built by a consortium of institutes and university 
departments from across Europe, Canada and the United States (NASA) under 
the leadership of SRON, Netherlands Institute for Space Research, Groningen, The 
Netherlands, and with major contributions from Germany, France and the US. Consortium 
members are: Canada: CSA, U. Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: 
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Observatorio Astronomico Nacional (IGN), Centro de Astrobiologia; Sweden: Chalmers 
University of Technology - MC2, RSS & GARD, Onsala Space Observatory, Swedish 
National Space Board, Stockholm University - Stockholm Observatory; Switzerland: 
ETH Zurich, FHNW; USA: CalTech, JPL, NHSC. Support for this work was provided 
by the Centre National de Recherche Spatiale (CNES), by the SCHISM project (grant 
ANR-09-BLAN-0231-01), by NASA through an award issued by JPL/Caltech, and by the 
Spanish MICINN (grants AYA2009-07304 and CSD2009-00038). This research has been 
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Table 1. LTE model parameters. 



Species 


Source size 


ex 


«LSR 


FWHM 


N 




(") 


(K) 


(km s -1 ) 


(km s -1 ) 


(cm -2 ) 


CH3OCH3 


33 


60 


-2.5 


5.5 


2 X 10 14 


CH3CN 


1 


150 


-1 


5.5 


7.4 x 10 14 


SO2 


1 


150 


-1 


5 


6.5 x 10 15 
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Table 2. HC1 column densities and abundances on the line of sight 

towards W31C. 



«->LSR 


7V(H I) a 


7V(H 2 ) b 


iV(H 35 Cl) 


iV(HCl) c 


A^HCl) / A^H d 


km s _1 


x 10 21 cm" 2 


X 10 21 cm- 2 


x 10 12 cm" 2 


x 10 12 cm" 2 


x 10- 10 


[15, 19] 


1.3 ± 0.35 


2.9 ± 0.3 


2.8 ± 0.51 


3.7 ± 0.67 


5.2 ± 1.1 


[19, 24] 


1.6 ± 0.56 


2.7 ± 0.3 


3.8 ± 0.69 


5.0 ± 0.91 


7.1 ± 1.5 


[25, 29] 


1.6 ± 0.40 


3.6 ± 0.3 


1.3 ± 0.28 


1.7 ± 0.37 


1.9 ± 0.5 


[30, 36] 


2.1 ± 0.78 


3.5 ± 0.3 


4.0 ± 0.74 


5.3 ± 0.98 


5.9 ± 1.2 


[36, 40] 


1.7 ± 0.20 


3.0 ± 0.3 


6.0 ± 1.02 


8.0 ± 1.3 


10.3 ± 1.9 


[15,40] 


8.4 ± 1.12 


15.6 ± 0.7 


17.8 ± 1.2 


23.6 ± 2.1 


6.0 ± 0.6 



a The H I column densities are obtained from I Fish et al.l 1 20031 ) absorption data with 
AT(H I) = 1.84xl0 20 (T Bpin /100K)/ r m dv [cm~ 2 /km s" 1 ], assuming a T spin - 100 K 



b The H2 column densities are derived from CH observations i Gerin et alj|2010 ). assuming 
iV(CH) = 3.5xl0~ 8 Af(H 2 ) jsheffer et al.lboosl) . 



c Assuming a solar 35 C1/ 37 C1 ratio of 3.1 
d 7V H =2A r (H 2 )+Af(H I) 
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Fig. 1. — Spectra of the J = 1 — transition of H 35 C1 (upper panel) and H 37 C1 (lower panel) 
towards W31C. The velocity scale in the upper axis is with respect to the frequency of the 
main (F = 5/2 — 3/2) hyperfme component of both isotopologues. The blue lines indicate 
the position of the HFS components. 
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Fig. 2. — Upper panel, XCLASS LTE models of the dimethyl ether (red), methyl cyanide 
(blue) and sulfur dioxide (green) emission lines blended into the H 37 C1 absorption line spec- 
trum, using the parameters given in Table 1. Additional lines (within the observations 
bandwidth) used in the model are shown in small windows at the top of the plot. The 
horizontal axis is in frequency units, with the corresponding velocity scale plotted on the 
top horizontal axis for consistency with Figure [U Lower panel, for illustration purposes 
only shows the resulting H 37 C1 spectrum (black) after subtracting the contamination from 
emission lines shown in red. The grey area highlights the emission line contamination-free 
area within the H 37 C1 absorption line spectrum. 
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Fig. 3.— LTE model of CH 3 CN (red), using CH 3 CN transition lines J = 29 - 28 (upper 
panel) at 533.123 GHz, J = 33 - 32 (middle panel) at 606.533 GHz and J = 34 - 33 (lower 
panel) at 624.878 GHz. The fit parameters are shown in Table 1. 



-23 - 



1.6 



0.6 



n 1 - 



n r 



n 1 1 - 




20 



40 

Velocity (km/s) 



60 



Fig. 4. — Spectra associated with each HFS components of the H 35 C1 line: F\ = 3/2 — 3/2, 
F = 2 - 2 (#reen), F x = 5/2 - 3/2, F = 3 - 2 (6/ U e); F Y = 1/2 - 3/2, F = 1 - 2 (red) 
(vertical offset introduced for clarity), the sum of the three components (magenta) and the 
observed spectrum (black) 
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Fig. 5. — Top panel, multiple Gaussian fit (red) to the hyperfine deconvolved spectrum of 
H 35 C1 J = — 1 (black) associated with foreground clouds on the line of sight towards 
W31C. Fo r comparison, we sho w the HFS deconvolved s pectra of Ho 35 Cl + (blue), H 2 37 C1 + 



(red) from 



Neufeld et al 



(120 121 ) and H I (magenta) from 



Fish et al. 



( 120031 ) 
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Fig. 6. — Left, abundance of the Cl-bearing molecules, chlorine atom or ion (X), normalized 
relative to the gas-phase chlorine abundance, /(X) = N(X.)/N(Cl), obtained with the 
Meudon PDR code, C1+ (magenta), CI (cyan), HC1 (blue), HC1+ (red), and H 2 C1+ (green) 
together with the H 2 fractions (black), /(X) = 2iV(H 2 )/(2iV(H 2 )+iV(H I)), averaged along 
the line of sight, as a function of total visual extinction A v tot . The initial density is 
tt-h = 316 cm -3 , and the radiation field intensity is 10. Right, comparison of the predicted 
fraction of gas-phase chlorine present in HC1, HC1 + , H 2 C1 + , CI and Cl + assuming a 
branching ratio of 90 % for chlorine and 10 % for hydrogen chloride (full lines) and the 
predicted fractions for a reduced branching ratio of 56 % for chlorine (dashed line). 



